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Abstract 

We perform a complete simulation of the process e + e~ — > t + t~vv, where v can be an 
electron, muon or tau neutrino, in the context of a general Higgs coupling to r leptons. 
We analyse various kinematical distributions and obtain the sensitivity regions in the 
parameter space that can be explored at a future e + e~ collider. In particular, inclusion 
of W boson fusion enhances the sensitivity significantly 
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1 Introduction 



The origin of fermion masses and mixings is one of most important issues in particle 
physics. Unfortunately, these parameters are inputs in the well-tested Standard Model 
(SM). All one can do is to measure them as accurately as possible and hope that in the 
future a more fundamental theory will actually be able to predict their values. 

Fermion masses are possibly related to the electroweak symmetry breaking mechanism, 
which is not known at the moment. In the simplest model, a scalar electroweak doublet 
in an ad hoc quartic potential is responsible for the symmetry breaking, leaving a scalar 
physical boson, the Higgs boson (J PC = ++ ), as a remnant. Favorite extensions of the 
SM, like the minimal supersymmetric standard model (MSSM), also predict the existence 
of a heavy pseudoscalar boson (J PC = h ), in addition to a light scalar boson. Another 
possibility is that the electroweak symmetry breaking is triggered by some new strong 
interactions, and in this case the lightest boson could be a pseudoscalar, like a pion [JEf. 
Therefore, it is very important to distinguish between these two cases. This can be 
achieved by studying the parity properties of the scalar particle, once it is found. The 
general framework for performing this study in the decay of the Higgs into fermion or 
gauge boson pairs was developed in and applied to Higgs production via Higgsstrahlung 
e+e- -> ZH in §. 

Possibly the most direct way to measure the CP properties of the Higgs boson is in 
the environment of a photon-photon collider, by setting the polarization of the photon in 
order to select different CP states ||. 

Details of the ZHH couplings can also be measured in a model independent way 
by studying the threshold behaviour and the angular distribution of the Higgsstrahlung 
process. In particular, the spin and parity assignements of the Higgs boson can be deter- 
mined p. 

At future hadron machines, like the CERN Large Hadron Collider (LHC), weak gauge 
boson fusion can provide an efficient way to measure the CP properties of the HWW 
coupling by analysing the azimuthal angle distribution of the two outgoing forward tagging 
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jets @. 

The coupling of the Higgs boson to top quarks can also be used as a measure of the CP 
properties of the higgs boson in the processes pp — > tth at the LHC ]7| and e + e~ — > ttH 
at a future Linear Collider ||. Prospects of measuring the parity of the Higgs boson have 
also been analysed in machines 0. 

In this paper we concentrate on the determination of the (pseudo)scalar-tau-tau cou- 
pling at a Linear Collider with a center-of-mass energy of ^/s = 500 GeV and accumulated 
luminosity of 1 ab" 1 , based on the Tesla proposal flQ |. We wil assume that this parti- 



cle has already been discovered at the Large Hadron Collider but a detailed study of its 
couplings is missing. 

Contrary to previous works []TT , 1^], where only Higgsstrahlung and Higgs radiation 



off a heavy fermion leg were considered, we take into account all relevant contributions 
to the process e + e~ — > t + t~vv, where v can be an electron, muon or tau neutrino. In 
particular, weak gauge boson fusion is the dominant contribution to Higgs production for 
M H < 180 GeV at yfe > 500 GeV. 

In extensions of the SM with extra scalars and pseudoscalars, the lightest spin-0 par- 



ticle can be an admixture of states without a definite parity ||12|| . Hence, we parametrize 
the general Htt coupling as: 

— (a + * 75 &) (1) 
v 

where v = 246 GeV and a = 1, b = in the Standard Model. We will present results 
considering a and b as independent parameters and also for the cases of fixed a = 1, free 
b and fixed 6 = 0, free a. We'll see that there is a region of insensivity around circles of 



\J o? + b 2 = 1 in the a — b plane. 

Without considering r decays, we are only sensitive to terms proportional to a , 
which comes from the interference with non-Higgs contributions, and a 2 and b 2 from pure 
Higgs contributions. Therefore, even at the r level we can search for deviations from the 
Standard Model prediction which could arise, for instance, in supersymmetric models. 

The analysis of r decays permits to study some P-odd correlations which isolates the 
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effects of the b term and hence can distinguish between a scalar and pseudoscalar (or a 
CP-violating mixture) nature of the Higgs-tau-tau coupling. This will be the subject of 
a forthcoming publication. 

2 Monte-Carlo simulation 

We generate all possible contributions to the processe e + e~ — > t + t~vv, where v can be 
either an electron, muon or tau neutrino. The following contributions are computed: 

• e + e~ — > T + T~v e i) e : there are 21 Feynman diagrams that contribute to this process 
(see Fig. [I]), of which only Diag. 10 (W-fusion) and Diag. 17 (Higgssthralung) are signal 
contributions; 

• e + e~ — > t+t^z/^z/^: there are 11 Feynman diagrams in this case (see Fig. [|), of 
which only Diag. 8 (Higgssthralung) contributes to the signal; 

• e + e~ — ► t + t~v t v t : there are 20 Feynman diagrams in this case (see Fig. ^), of 
which only Diag. 17 (Higgssthralung) contributes to the signal. 

The cross sections for these three contributions, with ran — 100 GeV for illustration, 
are respectively o Vt = 2.81 x 10~ 2 pb, a V/i = 1.25 x 10 -2 pb, and a Ve = 0.109 pb. As 
expected, a Ve is dominant due to the fusion diagram. 

In Fig. |] we show the dependence of the total cross section e + e~ — > t + t~vv, summed 
over the three neutrino species, on the parameters Aa = a — 1 and b. The dependence 
in b is symmetric since only terms proportional to b 2 contribute to the cross section. On 
the other hand, the dependence on Aa is asymmetric due to interference with standard 
model processes, which leads to the presence of linear terms in Aa. Our goal is to see 
how sensitive the experiments performed at the next generation of e + e~ colliders will be 
in the determination of these parameters. 

3 Analysis 
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3.1 Experimental Aspects 

To proceed with the analysis of the e + e~ — ► t + t^vv process one needs to have the final 
state topology as well as the r's momenta reconstructed. 

For the generated data sample of the e + e~ — > t + t^vv process we assume that tau- 
lepton pairs in the final state can be reconstructed in the experiment by using r decays 
into 7r and p. The r — > ttv and r — > pv decay channels are the two most important and 
together provide about 13% of the total branching ratio of r pairs. As it was shown in 
the LEP/SLC studies, other channels, such as r — > 7t~7t + tt ± i>, — > h°'k q 'k ± v 1 as well 
as leptonic tau decay modes, can be included increasing the experimental data sample to 
82% . 

There are several methods which can be used to reconstruct the final r pair topol- 



ogy fl3 |. Recently, a new method [14| was proposed, which allows one to have a high 
reconstruction efficiency in the — > p ± u T (u T ), p^ 1 — >• tt^tt decay chain. This method 
was used to study the Higgs boson's parity using r — > pv fll5| . 



The results of the LEP experiments [161 show that the r reconstruction efficiency 



varies from 51% to 80% . Based on the above results we adopt the moderate r-pair 
reconstruction efficiency e TT = 50% . 

Another important aspect is the assumption about the detector performance and pos- 
sible sources of the systematic uncertaintes. We include the anticipated systematic errors 
of 0.5% in the luminosity measurement, 1% in the acceptance determination, 1% in the 
branching ratios, and 1% in the background substraction, and assume the Gaussian na- 
ture of the sytematics. To place bounds on the Ht + t~ couplings, we use a standard 
X 2 -criterion to analyse the events. We found that the most strict bounds are achieved by 
dividing the distribution event samples into 10 bins. 

As for the possible contribution from background processes, like e + e" — > e + e~ ZZ — > 
e + e^T + T^vv (with e-pair lost), e + e~ — > vi)W + W~ — > uut + t^uu, e + e~ — > ZZZ — > 
t + t~vvvv etc., the cross-sections of these processes are either small, or they can be 
significantly supressed without much difficulties, as it was shown in the analogous study 
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for the process e + e — > vvH — > vvbb [I?]] 



3.2 Kinematical Distributions 

We analysed various distributions in order to determine which one is most sensitive to 
new physics: in the T correlation, defined by T = ^J-^ 3 p ei ■ (p Tl x p T2 ), the r momentum 
\p T \, the angle between r and the beam cos^ and the r pair invariant mass M T + T -. 
They are shown in Fig. |5] for the Standard Model, all plotted with 10 bins for illustration 
purposes. 

In Fig|| we illustrate the difference in the Higgs contributions in these distributions 
for the standard model (a = 1, b = 0) and (a = 0.5, b = 0.5). We see that the shapes are 
very similar, as expected, but the levels can be noticeable different. 

In order to find the level of sensitivity, we show in FigjT] the quantity ' — ^ or 
each bin of the distributions, where <j New is the cross-section with a and/or b different 
from their SM values. In this figure, we used a = 1, b = 0.5 as an example of new physics. 
A luminosity J Cdt = 1 ab _1 was used. We would like to stress again that including the 
contribution from W boson fusion enhances significantly the sensitivity and the result 
seems to be not so pessimistic as pointed out in [|ll|] . 

In order to determine the sensitivity regions in the a — b plane that can be probed 
in the next generation of linear colliders, we used the standard x 2 method where the 
experimental error Aa^ xp is given by: 



Aa^ xp = af M J 5 2 syst + 5 2 stat (2) 



where 



>stat 



af M e 



TT 



J Cdt 



(3) 



and 5 2 syst is the sum in quadrature of the systematic uncertainties mentioned above. 

In order to understand the relative contributions for each class of diagramas in Figs. 1- 
3, we show in Fig.[8]the individual contributions from electron neutrinos, muon neutrinos 
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and tau neutrinos channels. As expected, electron neutrinos dominate due to the weak 
gauge boson fusion contribution. 



4 Results 

In order to determine the sensitive region in the parameter space, we decided to use 
the cos^t- distribution, since the r direction is possibly the easiest observable to be 
reconstructed from data. We also found that using 10 bins in the analysis optimizes the 
results. 

We investigated the effects of the luminosity and the results are presented in Fig. ^|, 
where the sensitivity regions for J Cdt = 100 fb _1 , 1 ab^ 1 and 10 ab -1 are shown. We 
can easily see what can be gained from an increase of luminosity at TESLA. 



In Fig. [10], we show our final results for a Tesla-like environment, with luminosity 
of 1 ab -1 and center-of-mass energy of 500 GeV and for Mh = 100 GeV. We adopted 
a 50 % r-pair reconstruction efficiency and a systematic error of ~ 2% , in accordance 
with Tesla Design Report |10| . The allowed region for independent Aa and b parameters 



at 95% confidence level is the area between the solid circles. The horizontal band with 
long-dashed lines is the allowed region for the b parameter keeping a — 1. The two vertical 
bands with short-dashed lines are the the allowed region for the Aa parameter keeping 
b = (of course only the right band containing the SM point is physical). 
The bounds that can be obtained at 95% confidence level are: 

1. The case of two independent parameters: 

(0.70) 2 < (Aa + l) 2 + b 2 < (1.23) 2 . 

2. The case of b = and free Aa: 

-0.30 < Aa < 0.23 . 

3. The case of Aa = 1 and free 6: 

-0.7 < b < 0.7 . 
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These results can be easily scaled for moderate variations in the Higgs boson mass 
around 100 GeV by multiplying the bounds by a factor (M#/100 GeV) 2 . 



5 Conclusions 

We have performed a complete analysis of the sensitivity to new hrr couplings from the 
process e + e~ — > t + t^vv at the next generation of linear colliders. These new couplings 
are predicted by many extensions of the Standard Model. We showed that forthcoming ex- 
periments will be able to probe deviations of Htt coupling. For a Tesla-like environment, 
we are able to constrain the couplings in the region (0.70) 2 < (Aa + l) 2 + b 2 < (1.23) 2 
for independent a and b parameters. In order to disentangle the individual contributions 
of the parameters a and b, and to determine the CP nature of the Higgs boson, further 
analysis of r decay products is needed. This will be the subject of a separate publication. 
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Figure captions 



Figure 1: Feynman diagrams for the process e + e _ — > r + T~z/ e z/ e . 
Figure 2: Feynman diagrams for the process e + e~ — > t+t^v^v^. 
Figure 3: Feynman diagrams for the process e + e _ — > r + r _ z/ T z/ T . 

Figure 4: The total cross section e + e~ — > t + t~vu dependence on Aa (solid line) and 
b (dashed line). 

Figure 5: Distributions in T correlation, the r momentum \p T \, the angle between r 
and the beam, cos# er and the r pair invariant mass, M T + T -. 

Figure 6: Contribution of Higgs diagrams to the differential distributions in the 
e + e~ — > t + t~vu process for standard model (a = 1, b = 0, black dots) and a = 0.5, b = 0.5 
(crossed dots). 

Figure 7: Sensitivity functions aN ^~ e ^ M for each bin of the distributions, where as an 
example we took a = 1, b = 0.5 for new physics. 

Figure 8: Relative contributions to the sensitivity functions of electron neutrino (black 
dots), muon neutrino (rectangular crossed dots), and tau neutrino (crossed round dots) 
channels. 

Figure 9: Allowed regions for C = 100 fb" 1 (the domain inside the solid circle - there 
is no inner circle for this case), 1 ab _1 (domain confined by short-dashed contours), and 
10 ab _1 (long-dashed contours). 

Figure 10: The allowed region for independent Aa and b parameters at 95% confidence 
level is the area between the solid circles. The horizontal band with long-dashed lines is 
the allowed region for the b parameter keeping a = 1. The two vertical bands with 
short-dashed lines are the the allowed region for the Aa parameter keeping 6 = 
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